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Abstract 

We derive an expression for the determination of the apsidal angles that holds good for arbitrary central 
potentials. Then we discuss under what conditions the apsidal angles remain independent of the mechanical 
energy and angular momentum in the central force problem. As a consequence, an alternative and non- 
perturbative proof of Bertrand's theorem is obtained. 
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I. INTRODUCTION 



In 1873, J. BertrandyJ] published a short but important paper in which he proved that there 
are only two central fields for which all orbits radially bounded are closed, namely: The isotropic 
harmonic oscillator field, and the gravitational one. Because of this additional degeneracy - in 
the language of modem group theory we associated the unitary group f/(3) with the harmonic 
oscillator and the orthogonal group 0(4) with the gravitational potential - it is no wonder that 
the properties of those two fields have been under close scrutiny since Newton's times. Newton 
addresses to the isotropic harmonic oscillator in proposition X Book I of the Principia, and to the 
inverse- square law in proposition XI [2]. Newton shows that both fields give rise to elliptical orbits 
with the difference that in the first case the force is directed toward the geometrical center of the 
ellipse and in the second case the force is directed to one of the foci. Bertrand's proof is concise, 
elegant, and, contrary to what one may be led to think by a number of perturbative demonstrations 
that can be found in textbooks and papers on the subject, fully non-perturbative. As examples 



of perturbative demonstrations the reader can consult references 



lon-pertur 

m. 



We can also find in 



the literature demonstrations that resemble the spirit of Bertrand's original work as for example 
1^ 7]. All perturbative demonstrations and most of the non-perturbative ones, however, have a 
restrictive feature, to wit, they set a limit on the number of possibilities of the existence of central 
fields with the property mentioned above to a finite number and finally show explicitly that among 
the surviving possibilities only two, the newtonian and the isotropic harmonic oscillator, are really 
possible. 

Let us now outline briefly Bertrand's approach to the problem. In his paper, Bertrand initially 
proves by taking into consideration the equal radii limit that a central force /(r) acting on a point- 
like body able of generating radially bounded orbits must necessarily be of the form 

/(r) = /cr(i/^''-3), 

where r is the radial distance to the center of force, k is a constant and p a rational number. Next, 
making use of this particular form of the law of force and considering also an additional limiting 
condition, Bertrand finally shows that only for p = 1 and p = 1/2, which correspond to Newton's 
gravitational law of force 
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and to the isotropic harmonic oscillator law of force 



/ (r) = -Kr, 



respectively, we can have orbits with the properties stated in the theorem. Moreover, Bertrand can 
also prove that only for these laws of force all bounded orbits are^closed. For more details the 
reader is referred to the original paper or to its English translation [1]. 

In the present paper, we consider the issue from an alternative point of view. This alternative 
approach stems from a different kind of question that we can ask ourselves. For radially bounded 
orbits there are two extreme values for the possible radii, namely, a maximum and a minimum one. 
These extreme values were named by Newton himself superior apse and inferior apse, respectively, 
or briefly, apses. The angular displacement between these two successive points or apses defines 
the apsidal angle A9a. In a central field of force, the apsidal angle depends on the total mechanical 
energy E and the magnitude of the angular momentum i . The question we can formulate is the 
following: For what potentials the apsidal angle has the same value for all orbits - that is for 
arbitrary energy and angular momentum - and what is its value? The answer as we will show in 
the next sections is that the apsidal angle is constant for only two potentials: The Newtonian and 
the isotropic harmonic oscillator one, and consequently as a corollary we have an alternative and 
non-perturbative proof of Bertrand's theorem. 

II. THE APSIDAL ANGLE 

In a central field, in which the magnitude of the force f (r) depends only on the distance r = ||r|| 
to the center of force, we can introduce the potential function V{r) with the property /(r) = 
—W{r) such that the total mechanical energy of a particle with mass m orbiting in this field is 
a constant of motion. Thus we write 



where v is the velocity of the particle. Moreover, the conservation of the angular momentum I of 
the particle in a central field constrains its motion to a fixed plane and allows for the introduction 
of an effective potential defined by 



(1) 



U{r-) = V{r) + 



2mr 



,2' 



(2) 
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with the help of which it is possible to reduce the motion to an equivalent unidimensional problem. 
Let us define for convenience the variable z = Xjr and the functions f (z) = ViXj z) and = 
IJ{\I z) in such a way that the effective potential Q is now written 

u{z) = v{z) + —z\ (3) 

Let us also assume that the effective potential ([3]) has a minimum at z^ = I/tq, such that the 
particle can move in a circular orbit with radius tq. This particular value of z is determined by the 
condition 

u\zo) = v'izo) + —zo = 0, (4) 
m 

from which we obtain 

/2 = -mv'{zo)/zo). (5) 
Taking this result into the effective potential ([3]) we can write 

u{z) = v{z) - — z . (6) 

ZZq 

In this way we see that the first derivative of the effective potential at the point z = zq is zero for 
each and every potential function v{z). This means that the vanishing of this derivative does not 
impose restrictions of any kind on the potential function v{z). The second derivative at the point 

Z = ZqIS 

u"{zo) = v"izo) - (7) 
Zq 

If u"{zo) is zero for an arbitrary point zq, we see that the potential must be of the form v{z) = 
a + bz"^, and in this case the force must be of the form n/r^. Newton in the Proposition IX of Book 
I of the Principia shows that in this case the orbit is an equiangular spiral and therefore without 
apsidal points. Hence, in what follows we will suppose that the second derivative of the effective 
potential is not null. 

If now we assume that u and z are complex variables and perform an analytical continuation of 



the function u(z) we can determine the inverse function z = z(u) 
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of the Biirmann-Lagrange series for a multivalued inverse function [llOh . With this purpose in mind 



applying the generalization 



we choose the point mq = u{zo) as the point around which we will perform the expansion of this 
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function. The result is 



Z = Zo 



n=l 



d 



where 



x(C) 



- X" (C) 



1 ) 
I 2 



(8) 



(9) 



and Mo = u{zo) is a first order algebraic branching point of the function z{u). It follows that we 
can define two real inverse functions of the function u = u(z); the first one is 



En 

n=l 

which holds for z > zq (orr < tq). The second one is 

oo 

Z2 = Zo + '^ (-l)"C„(u - Uo) 2 , 



(10) 



(11) 



n=l 



which holds for z < Zq (or r > tq), and where we have written for convenience 

1 d^-'x^iC) 



Cri 



n-l 



(12) 



It is also convenient to write the function ^ in the form 

1 



x(C) 



fv{C) - v{zo) - v'{zq) (C - Zq) v'{zo) 



(13) 



(C - z,f 2zo 

which can be obtained by making use of ([6]). For radially bounded orbits the two extreme values 
for the radii, namely, r^ax and r^^^, maximum and minimum (or z^^^ and Zmax)^ respectively are 
determined by the condition = 0, or Za = 0. The particle oscillates indefinitely between rmax 
and rjnin- In terms of the effective potential radially closed orbits are characterized by extreme 
points that satisfy the condition E = U{ra), or E = u{za). For convenience we take the direc- 
tion defined by the arbitrary vector tq as the reference for the measure of angular displacements. 
Hence, the angular displacement between two successive apses, that is, the apsidal angle A 9a, 
can be written in the form 



(14) 



where A 9i is the angular displacement from the point to the point r^i^ and A 6*2 is the angular 
displacement to the apsidal point rmax- 



The angular displacement in a central field of force can be easily determined from the conser- 
vation laws of the mechanical energy and the angular momentum. Hence we can write 

i dr 



A9i 



ro 



mr^ / 2 



ra 



{E-U) 



and 



dr 



mr^ / 2 



-{E-U) 



m 



Making use of the inverse functions (flOl ) and (fTTl) we have 



dzi dU 



^ Juo dU ^{E - U) 



^nCn 



and 



Mo 



dz2 dU 



dU, 



2m Juo dU y/(E - U) 



The integrals in equations (fTTl) and (fTSi) can be evaluated with the help of 



/ (x — a)^ ^ {b ~ xY ^ dx = {b — a)^^" ^ B (A, z/) 

J a 



where 



(15) 



(16) 



(17) 



(18) 



V{\ + z/) 



which holds for 6 > a, 3ft /i > and 3ft z/ > 0, see [11], Formula 3.196.3. Therefore we have 



2(2A;)!! 



k=0 

oo 



(E-U,f^\ 



k=0 



and 



A9n 



in 



'i2k + l)\ 
^ (2k -IV^ 



2V2m 



k=0 



(2A;+ 1)!! 



(19) 



{E-Uo 



(20) 
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Notice that the first term in the first summation is the only one that does not depend on the energy. 
Notice also that the constants Cn do not depend on the energy as well. After adding equations (fT9l) 
and (|20l ) we obtain for the apsidal angle the expression 

^»" = ;|;jEP* + i)c^-«^7S)if(^ (21) 



Notice that this series does not depend on the even coefficients. Equations (fT9l) . (1201) . and (|2T]) 
are our main results and in the next section we will use them to answer the question we posed at 
the beginning of this work, to wit: For what potentials the apsidal angle has the same value for 
all orbits, i.e.: for arbitrary values of the energy and angular momentum, and what value does it 
assumes. 

III. THE NEWTONIAN POTENTIAL 

Let us begin by looking for potentials that keep the partial angles A 9i and A 62 separately 
constant. If this happens to be so, the apsidal angle A 9a will also be constant for any value of the 
energy and the angular momentum. In order to accomplish that it is necessary that for n ^ 2 all 
coefficients C„ in equations (17) and (18) be equal to zero. After evaluating the coefficient C2 and 
setting it equal to zero we obtain v"'{zo) = 0, and taking into account that zq is an arbitrary point 
it follows that 

v"'{z) = 0. (22) 



The general solution of (|22|) is 



v{z) = az^ + bz + c. (23) 



The integration constant c is an additive term to the potential and can be discarded without loss of 
generality. The second derivative of the effective potential at the stationary point is u" = —h/zQ. 
It follows that Zq is effectively a minimum only if the constant h is negative and this corresponds to 
an attractive field. For a potential function of the form given by (l23l) we can evaluate the function 
x(C) given by (fT3l) to obtain 

X(C) = -1^- (24) 




2z, 

Hence the only non-zero coefficient is Ci, because the function x(C)' for this potential, does not 
depend on ^ and therefore A6'i = A 6*2 = constant. Evaluating the first coefficient for this potential 
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we obtain 

Ae, = Ae, = l^i + ^-^. (25) 

For a fixed value of tlie angular momentum any potential of the form given by (l23l) generates 
bounded orbits with equal and constant angles A 9i and A 62 and thus the apsidal angle A 9a as 
defined by (fT4l) will also be constant. Moreover, if only the constant a vanishes, that is, if the 
potential is of the form v{z) = bz, it is possible to obtain a potential for which the apsidal angle 
is also independent of the angular momentum. This potential corresponds to the newtonian one, 
the apsidal angle is vr and the orbits are closed. The vector tq is perpendicular to the vectors Fmax 
and Fmin, which in their turn are antiparallel vectors with respect to each other defining in this way 
only one symmetry axis of the trajectory of the particle. 



IV. THE ISOTROPIC HARMONIC OSCILLATOR POTENTIAL 



Let us now look for potentials that keep the apsidal angle A6a constant. In order to accomplish 
this it is necessary that all the odd coefficients C„ starting from the third one be zero. As before we 
begin by explicitly calculating the lowest order coefficient, namely C3. Evaluating this coefficient, 
setting it equal to zero, and considering the arbitrariness of zq we obtain the following differential 
equation 



5 ///2/ N 



v"{z) 



v'{z) 



v^'^'Hz) 



0. 



Writing the function v{z) as 



v{z) = - 



j ^^^^ ^\ j 



we can recast (|26l) into the form 



1 d 



[z^{z)] - 



3zv{z)dz'~^'~'' lilz^iz)] 
Equation (|28l ) can be immediately integrated and its general solution is 



(26) 



(27) 



(28) 



1 



{Az' + B) 



-3/2 



(29) 



where A and B are arbitrary constants. Notice now that the newtonian potential analyzed before, 
which is a solution of (l22l) . is also necessarily solution of (|26l ) as can be immediately verified. Such 
a solution can be obtained from (|29l) by setting A = 0. Therefore, without any loss of generality 
we suppose A 7^ and for convenience recast (|29l) into the form 
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<P{z) = ^{z' + b')-'/\ (30) 

where k ^ and b are new constants. We write 6^ in order to assure that the function 0(2;), 
and consequently the corresponding potential function, does not become ill-defined in the region 
around z = 0. Making use of (|30l) and (ITTI) we can determine the potential function v(z). For 
6 = we obtain 

k 

^iz) = — + cz^ + d, (31) 



2z 



and for b ^ 



v{z) = + cz^ + d. (32) 

The constant c? is a simple shift of the zero of the potential and does not influence the law of force. 
Considering first the potential given by (|3T1) and evaluating the second derivative of the effective 
potential at the stationary point we obtain 

u"{zo) = -i. (33) 
^0 

It follows that we will have a minimum only if A; > 0. Making use of (fT3]) we have 



Performing the analytical continuation of (1341) and making use of the Cauchy formula we obtain 
from([l2l) 

/ \ — 

^^ = -1?) T^^o" / ,,n dC. (35) 



n \k) 27,1 ' Jc^^ (C^-z, 



^0. 



The integral in (|35] ) can be easily evaluated by making the transformation r] = ('^ — Zq in the 
neighborhood of Zq. It follows then 



The residue of (l36l) is obviously null if n is an odd integer greater than one. Therefore, we can 
be sure that all potentials given by (|3T| ) generate orbits for which the apsidal angles depend only 
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on the angular momentum. The apsidal angle for this class of potentials follows from (121]) and is 
given by 



In order to have an apsidal angle not dependent on the angular momentum it is necessary that the 
constant c be set equal to zero. This will lead us to potentials of the form 



^i^) = —2^ (38) 



k 

which corresponds to the isotropic harmonic oscillator. In this case (|37l) yields 



A^a = |, (39) 



and, therefore, the orbit is closed. 



V. THE ADDITIONAL THIRD POTENTIAL 



Finally let us consider the potential given by (|32l) for which the effective potential reads 



u \z) = — , (40) 

^ ' Ab^zi^/4T¥ 

The second derivative of u{z) at zq is easily evaluated and the result is 

u"{zo) = , " (41) 

We can see that only for a > we will have a minimum and in what follows will show that this 
must indeed be the case here. Making use of (fT3l) we obtain after some algebraic manipulations 



1 



(62 + 2^2) e + 2^4 + 62^oVC + + 



x(C)-V^VV^o + ^^^o ^ . (42) 

Next we consider the analytical continuation of this function. The complex function given by (l42l) 
has algebraic branching points at z = ±zo and z = ±ib and is defined on a four-sheeted Riemann 
surface. Hence we choose a sheet on this Riemann surface such that the functions x(C) does not 
have branch points at ^ = ±Zo and with a branch cut along the imaginary axis from z = ib to 
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z = ~ih. Notice that in this case the numerator of (l42l) does not have zeroes. Making use of 
Cauchy formula we can write 



dz 



n-l 



{n-l)\ 



z=zo 



2711 



^" (C) dC, 



(43) 



where we have defined 



' = V 2* ^ —I 

and is a small closed path encircling the point Zq. The integrand in equation (|43l) has another 
pole at the point z = —zq. Then, choosing the path indicated in Figure ^ the following results 
can be easily obtained 



lim /?9"(C)ciC = 0, 

i?— >oo J 



(45) 



if n > 1; 



hm /^"(CK = 0, 



and 



hm lr{C)dC = 0. 



(46) 



(47) 



Finally, the integral over the finite parts of the real axis and the imaginary axis, see Figure 
cancel out. Then we have 



^^{C)dC+ j r(C)rfC = 0, (48) 
(=zo C=-2o 

for n > 1. We can easily prove by making the substitution ( —( in the second integral above 
that 

+ J ^"(C)c^C = 0, (49) 

and therefore we conclude that 



^n-l^n (Z) 



(50) 



z=zo 
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if n is an odd integer greater than one. The last step is the evaluation of the apsidal angle for this 
case. Making use of (01) and (140)) we have 



m 



mv' {zq) m 

Zq Zq 

a {2zl + b) 



1 a {Azl + 2hzo 
'2 h^^^Thzi 



2c 



h'^ZQ^Jzl + h 
Taking this result into (|2TI) we obtain 



2c 



(51) 




62 



zl + b 



(52) 



The potentials given by (|32l) generate orbits for which the apsidal angles depend only on the 
angular momentum. In this third case it is not possible to find constants h and c in such a way that 
those angles do not depend also on the angular momentum. 



VI. FINAL REMARKS 

In this work we have introduced a novel method for the determination of the apsidal angle in an 
arbitrary central field of force. For this reason equation (|2T1) should be considered as the main result 
of this paper. From this equation we have discussed the conditions under which the apsidal angles 
remain independent of the energy and angular momentum and consequently lead to bounded and 
closed orbits, and we have found that only the Newtonian and the isotropic oscillator potentials 
present such behavior. Moreover, from equation (|2TI) we have also calculated explicitly the value 
of their associated apsidal angles and have found the well-known results 7r/2 and tt, respectively. 
We can see also that these two special fields are strongly determined by the dependence of the 
centrifugal potential on l/r^. As a consequence, we have re-obtained Bertrand's theorem in an 
alternative but non-perturbative way. After more than one hundred years since its publication 
Bertrand's theorem still fascinates us uM- Any serious student of classical mechanics sooner or 
later will become at least intrigued by it. No wonder that alternative proofs some of them very 
interesting - see for instance the phase space approach in lll3n or the hamiltonian one in [|l4ll - 
were offered in the past and nowadays. We have also obtained a third potential which we can lay 
aside because it does not meet the conditions that lead to bounded closed orbits. 
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The analytical function techniques applied to the problem of finding the only central fields that 
allow an entire class of bounded, closed orbits with a minimum number of restrictions lead in a 
concise, straightforward way directly to the two allowed fields. It is also worth to mention that, 
as far as the present authors are aware of is one of the few examples of the usefulness of the 
Biirmann-Lagrange theorem in its most general formulation. Equation (|2T]) seems promising and 
its applications go beyond the rederivation of Bertrand's results. In fact, it can be taken as the 
starting point in the discussion of problems related to precession phenomena. Work by the present 
authors in this direction is in progress and results will be published elsewhere. 
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FIG. 1 : General form of the effective potential energy. 
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FIG. 2: Contour 
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